Effective charges and statistical signatures in the noise of normal 
metal— superconductor junctions at arbitrary bias 
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Shot noise is studied in a single normal metal-superconductor (N-S) junction at finite frequency, 
and for branched N-S junctions at zero frequency. The noise spectral density displays a singularity 
at the Josephson frequency (ui — 2eV/h) when the applied bias is smaller than gap of the super- 
conductor. Yet, in the limit eV S> A, quasiparticle contributions yield a singularity at u) = eV/h 
analogous to that of a normal metal. The crossover between these two regimes shows new structures 
in the noise characteristic, pointing out the failure of the effective charge model. As an alternative 
to a finite frequency measurement, if a sinusoidal external field is superposed to the constant bias 
(non stationary Aharonov-Bohm effect), the second derivative of the zero frequency noise with re- 
spect to the voltage exhibits peaks when the frequency of the perturbation is commensurate with 
the Josephson frequency. Finally, the statistical aspects of noise are studied with an analog of the 
Hanbury-Brown and Twiss experiment for fermions: a superconductor connected to two normal 
leads. Noise correlations are found to be either negative (fermionic) or positive (bosonic), due to 
the presence of evanescent Cooper pairs in the normal side of the junction, in the latter case. 
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I. INTRODUCTION 



Recent developments in condensed matter physics em- 
phasize the importance of shot noise in mesoscopic con- 
ductors. Noise contains more information than the 
conductance: in the Poisson regime for instance, zero- 
frequency noise is proportional to the current and to 
the effective charge. A convincing application is the 
direct measurement of the fractional charge in a quan- 
tum Hall liquido. Another important issue of noise con- 
cerns the statistics of the (quasi-) particle^ It is now 
established that Hanbury-Brown and Twisso type corre- 
lation experiments yield a different sign for bosons and 
for fermionsaij. 

This paper focuses on both issues in normal metal- 
superconductor (NS) junctions. Existing results on junc- 
tions between two normal metals (NN) are summarized 
below. Finite frequency shot noise has a Angularity at 
Tllo = eV (where V is^the applied bias)Q, which was 
detected experimentally!] . Another phenomenon, analo- 
gous to a finite frequency measurement, called the "Non- 
Stationary Aharonov-Bohm effect"EI uses a local alter- 
nating field superposed to the bias voltage. Steps in 
the noise|-iderivative with— respect to the DC bias were 
predicted^ and measuredliiJ. The height of these steps is 
non-monotonic with the amplitude of the harmonic per- 
turbation. 

The statistical signatures of noise correlations can be 
illustrated in a three terminal device where current fluc- 
tuations in the twa,ceceiving leads are expected to be 
fully anti-correlatecHu, a direct consequence of the Pauli 
principle. Experiments were performed in the integer 



quantum Hall regimeQ and in a two dimensional electron 
gasQ, with a splitter used to partition an incident beam 
of electrons into reflected and transmitted beams. The 
measurement of the correlations between the reflected 
and the transmitted beams reported a negative value, 
confirming the theoretical predictions. 

In contrast to normal junctions, transport in NS junc- 
tions involves Andreev reflectionEa: an incoming electron 
is reflected as a hole at the boundary. The remaining 
charge 2e is absorbed by th&,superconductor via a Cooper 
pair. In a previous paperl!3, finite frequency noise was 
computed in the Andreev regime only (eV <C A). The 
noise characteristics then present a singularity at the fre- 
quency oj — 2eV/h, a signature of an effective charge 2e 
corresponding to that of a Cooper pair. This issue is not 
surprising, because in the Andreev regime, some phys- 
ical quantities can, in principle, be extracted from the 
normal metal results by replacing e by 2e. In particu- 
lar, at zero-frequency the doubling of shot noise and the 
crossoverJxam thermal noise to excess-jipise have been 
predicteal3il3 and recently measurecOEZl. However, for 
thermal noise (at ksT -C A), this naive substitution 
would imply a violation of the fluctuation dissipation the- 
orem. 

Failures of this effective charge model also occur when 
the bias is comparable to the gap or when it is increased 
beyond it. Indeed, when eV > A, the effective charge 
transfer picture breaks down, as quasi-particles above 
the gap bear a single electron charge. It will be shown 
that at such voltages, the cusps/singularities which ap- 
pear in the finite frequency noise neither correspond to a 
charge e nor to a charge 2e. Thus, a careful analysis and 
an explicit calculation are especially needed to describe 
the crossover from sub-gap to above gap regime. More- 
over, with the advent of superconducting samples with 
a "small" gap, this also allows to consider the limit of 
large biases (eV ^> A), where single quasiparticle trans- 
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mission overrides Andreev reflection. In this situation, 
one recovers the results obtained for the NN junctions, 
with a charge transfer e. 

At the same time, an NS junction which contains a 
splitter on the normal side constitutes a Hanbury-Brown 
and Twiss tepe experiment where statistical effects can 
be detectedta. In the Andreev regime, the noise corre- 
lations could possibly be positive (bosonic) due to the 
presence of Cooper pairs on the normal side (proximity 
effect), whereas a quasi-particle dominated regime should 
favor negative (fcrmionic) correlations. The purpose is 
here to investigate the effective charges and statistical 
tendencies which show up in the noise characteristics of 
an NS junction in an united fashion. Results for the sub- 
gap (Andreev) regime will be recalled, and confronted to 
novel results for above gap transport. 

The paper is organized as follows. Noise correlations 
in a multi-terminal device coupled with a superconduc- 
tor are computed at finite frequencies (section |l|). In 



section [II a single N-S junction is stu died for several 
cases: ajm the Andreev regime (section |IIIB ), the noise 
spectral density presents a singularity at the Josephson 
frequency; b) wh en the applied bias is increased beyond 
the gap (section III C ), additional singularities appear; 
c) if a sinusoidal external field is ad ded (non-stationary 
Aharonov-Bohm effect, section III D ), the second deriva- 
tive of noise with respect to the bias presents peaks when 
the frequency of the perturbation is commensurate with 
the Josephson frequency. The last section ( IV ) deals with 
the fermionic Hanbury-Brown and Twiss experimental 
proposal with a superconductor, showing that noise cor- 
relations may be either negative or positive. 



II. CURRENT, NOISE AND NOISE 
CORRELATIONS 



A. Assumptions 

The superconductor is connected to an arbitrary num- 
ber of normal leads (multi-terminal device), and the de- 
vice is supposed to be small enough that all scattering 
processes are elastic. Note that the case-af two supercon- 



ductors is not addressed here (see Ref.td). Calculations 
are restricted to the one-channel case, but a generaliza- 
tion to a multi-channel system is straightforward. Thus, 
transport is dominated by the scattering properties of the 
N-S junctionnJCJ. The two scattering processes at play 
are then Andreev reflection and normal reflection of elec- 
trons and holes. The superconductor is maintained at a 
constant chemical potential fig, and each normal termi- 
nal is fixed at the same potential hn. For convenience, 
all energies are measured with respect to fig: so that the 
applied bias reads /txjv — Ms = This bias is chosen to 
be positive throughout the paper. 

Shot noise in a given lead, or alternatively noise cor- 
relations between two (normal) terminals are defined as 



the Fourier transform of the current-current correlation 
function: 



Sij(u>) = lim 



f 



T^+oo T 



T/2 



-T/2 



dt I dt'e lult ' 



((/iW/^t + f))-^)^)) ■ (1) 



When i — j, Su(lo) corresponds to the noise in the ter- 
minal i, whereas if i and j are different, Sij(O) represents 
the zero frequency noise correlations between lead i and 
lead j. ( ) designates the thermodynamical average in 
the grand-canonical ensemble. 



B. Bogolubov-de Gennes equations 

The Bogolubov-de GenneslEI (BdG) approach to inho- 
mogeneous superconductivity is the adapted formalism 
to treat electrons and holes on the same footing. Per- 
forming the Bogolubov transformation (which must di- 
agonalize the effective BdG Hamiltonian) , and going to 
an energy representation, the annihilation operator of a 
particle with spin a (a = ±1) at the position x in the 
terminal i ipi i(T (x) can be written as: 
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'IE I Uiip{x) c ]P<J {E) 



'hvi(E) 



—a 



«<j/g(s) 



-j -a- 



(E) 



(2) 



State Uijf3 (vijp) corresponds to the wave function of a 
electron (a hole) scattered in terminal i, due to a quasi- 
particle of type (3 (electron or hole, /3 = e,h) which was 
incoming from lead j. Operators c(E) and c + (E) satisfy 
standard anticommutation relations. v J e (E) — fik J e (E) is 
the velocity in the lead j. The Bogolubov-de Gennes 
equations may be written as: 



Eu ij!3 {x) 



Us + V(x) Uijp(x) 



Evijp(x) = - I - 



Ira dx 2 

+A*(x)u lj p{x) , 



(3) 

and describe the evolution of particle states Uijp and 
Vijp. The pair potential A(x) should be calculated self- 
consistcntly, but for simplicity corresponds here to the 
superconducting gap in the bulk superconductor (x > 0), 
and gives zero in the normal terminals (x < 0). 
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C. States in a normal terminal 

In order to calculate noise with Eq. ([!]), the states 
which appear in the Bogolubov transformation (|^) are 
specified using the scattering (5) matrix describing the 
junction. In a normal, ideal lead A(x) — and V(x) = 0, 
the Bogolubov-de Gennes equations (|J) reduces to a 
Schrodinger equation for electrons, and to its time re- 
versed analog for holes. Solutions of the form e = x 
for electrons, and e ^ x for holes are chosen, where 
fcf = ^/2m {p s + E)/Ti and fcf = ^J2m - E)/fi are 
the wave vectors of electrons and holes. Electrons and 
holes in lead i which originate from a particle of type f3 
(f3 = e, h) in terminal j and scattered into i are described 
by: 



ijef3 




(4) 



v ijP (x) = SijS^pe lk * x + s ijhfi \l T^e^'* . (5) 

Note the opposite sign for momenta of electrons and 
holes. Xi, the position in terminal i is specified as x in 
Eqs. (||) and (||) for simplicity of notation. Sij a p is the 
scattering matrix element expressing the amplitude of an 
outgoing particle a in lead i due to an incident particle 
of type (5 in lead j. 



D. Current operator and average current 



The current operator in lead i is defined by: 



dip iltr (x) 



dx 



dx 



(6) 



Substituting ip by Eq. (|2|), the current operator becomes: 
eh 



Ivnivp 2irh , 



i r-t-oo /"too 



^u irn d x Ui n d x u irn Ui n j< ln rr < ,t 



+ {v im d x v. 



i7n yj x u in ^x^im^in 1 Cm — a — a 



(7) 



where the sums over j and (3 have been replaced by 
a single sum over index m. Expressions with index 
m (n) have an energy dependence E\ {E%). Solv- 
ing the Bogolubov-de Gennes equations in the su- 
perconductor leads to find the waves vectors kf h = 



1/2 



y/2m(fi s ± (E 2 - A 2 ) 1/2 J /% for electron-like (+) 

and hole-like (-) quasi-particles. The chemical potential 
of the superconductor [1$ is large compared to all the en- 
ergy scales involved in the system. Thus, the assumption 
kf = k^ = A;f = kfl = kp is made, which explains the 
presence of the Fermi velocity vp in Eq. (Q). 

The calculation of the average current involves the av- 
erage of creation and annihilation operators products 
such as (c+ rJ (E 1 )c niJ (E 2 )) = f m {Ex)5 mn S{Ei - E 2 ). 
f m (E) = Jfd(E =f eV) for electrons and holes on the 
normal side, with fpo(E) the Fcrmi-Dirac distribution. 
fm(E) = Jfd(E) in the superconductor. The average 
current is then: 



(h(x)) = 



2nmivF 



r+oo 

-(v im d x v* m - d x Vi m v* m ^ (1 - f m ) . (8) 



E. Noise and noise correlations 



Since a single current operator is composed of products 
of two creators or annihilators, (Ii(x,t)Ij(x,t + t')) is a 
sum of average values of four creation or annihilation op- 
erators. These average values are expressed as a function 
of the Fermi-Dirac distributions using Wick's theorem. 
The calculation of noise is now performed using Eq. (Q). 
It is convenient to define the following matrix elements: 



Mm j 



n {E,E',t) = u jn (E',t)d x u* im (E,t) 

-u* im (E,t)d x u jn (E',t) , (9) 
B imjn (E, E', t) = v* n (E', t)d x v m (E, t) 

—Vim(E,t)d x Vj n (E',t) , (10) 

Ci m jn(E, E ,t) — Uj n (E' ,t)d x V im (E,t) 

(E,t)d x u jn (E',t) . (11) 

Calculating all the average values and the difference 
(Ii(t)Ij(t + t')) — one obtains the noise or the 

noise correlations: 



2*2 



1 



lim 



2m 2 v F (2<Khy T J_ T/2 



+T/2 

dt 



dt'e iu>t ' 



'- + CJO 

: / dE 
'o Jo 



dE'J2< 



fm(E)(l-f n (E'))e^ E '-^^ h 

x A imjn (E,E',t)A* mjn (E,E',t + t') 
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+ B imjn( E ' E ', t)B lmjn (E, E', t + t') 

+Aimj n (E, E , t)B im .j n (E, E , t + t) 

+ B imjn( E i E ', *)^*mjn(- E ) E> , t + 

+.f m (E)f n {E')e-^ E + E 'V' h C* mjn (E, E') 
{E',E,t + t r ) 

+Ci m j n (E, E' , t + t')j 

+ (1 - f m (E))(l - f n (E'))e^ E + E '^/ n 
x (c ]n i m (E' , E, t) + Ci m j n (E, E 1 , i)^ 



(12) 



+e(-E - hu)(i - f m {-E - huj)){\ - f n (E)) 

xC* min (-E-hu,E) 
x(c inim (E,-E -Tiw) 



+Ci m i n ( — E — hui, E) 



(14) 



The zero-frequency limit, .for noise correlations between 
terminals i and j yieldsE§Eai3E£l: 



e 2 h 2 1 



m 7 2m 2 v| 2 



1 / >+oc 

— / d£^/ ro (£)(l -/„(£)) 

«/0 ™ „ 



x I A imjn (E, E) + B* imin {E,E)\ . (15) 



So far as the process is stationary (for example if no 
time dependent external field is applied), matrix elements 
Aimjn, E imjn et Ci m j n are time independent. As a con- 
sequence, the integration over t simplifies, and the inte- 
gration over t' gives 8 functions with energy. As a result, 
terms proportional to (1 — / m )(l — /„) cancels because 
one is dealing with quasi-particles with a positive energy. 
This yields: 



Sij(u>) 



e 2 h 2 
2m 2 



vl2irhj 

* u m.n \ 



Q(E + hu)f m {E + hu)(l - ME)) 
x | Ai m jn {E + hu,E) + B* mjn {E + Tiuj,E)\ 
+e{hto - E)f m (hu - E)f n {E) 
xC; mjn (huj-E,E,t) 



hui - E) 



+Cimjn{TT'U — E, E) 
-Q(-E - M(l -fm{- E - M)(l -fn( E )) 

xC: min (-E-hu,E) 



X (C)ni m ( E , —E — Huj) 



Cimjn{ — E — TlU, E) 



(13) 



after integration over E' . This expression corresponds to 
the noise if indices i and j are the same: 



Su(uj) 



3 2 h 2 



1 



dE 



2m 2 Vp 2-kTi j 
Q(E + hu)f m (E + M(l - fn(E)) 
x \A lmm (E + hu, E) + B* mm {E + hu>, E) \ 

+e{hLo - E)f m (hLo - E)f n (E) 

xC; mm (hio-E,E,t) 

X (Ci n i m (E, hut — E) + CiminQlU) — E, E) 



III. SINGLE N-S JUNCTION 

A. General expression for noise 

A single N-S junction with arbitrary transparency 
(with a tunneling barrier at the interface) is considered 
first for a small applied bias (Andreev regime eV <C A) 
and then for biases larger than the gap. Analytical ex- 
pressions are obtained in the first case, while numerical 
results will be presented in the latter regime. For simplic- 
ity, the temperature is chosen to be much smaller than 
the gap in both cases. In the the Andreev regime, the 
integrals over energy in Eq. ( |l4| ) can be performed, re- 
sulting in three distinct contributions combining different 
products of Fermi-Dirac distributions. It is interesting to 
note that although the current operator of Eq. (Q) itself 
cannot couple states which differ by two quasi-particles, 
its fluctuations give a contribution which is proportional 
to fnfm in this finite frequency calculation. Rewriting 
down A mni B mn and C mn as a function of the S ma- 
trix elements, three distinct expressions of the noise are 
found. If hus < eV: 



o 2 f reV-hu 



\s NNee (E + huj)\ 2 (l - \s NNee (E)\ 2 ) 

+ \sNNhe( E + ftu)\ 2 (f " \s N Nhe(E)\' 
+ s *NNeei E + ^) s N N ee( E ) 

X S N Nhe(E + huj) s* NNhe (E) 
+s NNee (E + hu>) s* NNee (E) 

x s *NNhe( E + hu) S NNhe (E) 



dE 



\s NNee (E)\ 2 \s NNhe (huj - E)\' 



+SjyNee(fou — E) SNNhe( E ) 

X s *NNee( E ) S* NNhe (huj - E) 



(16) 
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If eV < huj < 2eV, one obtains: 
dE 



2e 2 reV 

= X , 



s NNee (E)\ 2 \s NNhe (huj - E)\ 2 

+ SNN ee (?LUJ — E) SNNhe{E) 
*S*NNee( E ) s *NNhe(Tw ~ E ) ■ ( 17 ) 



If hu> > 2eV noise vanishes. Eq. fll7| ) contains no spa- 
tial dependence as a consequence of the approximation 
^5 3> Tiw, mentioned in Sec. ( |II C| ). 



argument for an oscillatory time dependence can be re- 
peated, since the incoming electron wave function oscil- 
lates as tp e ~ exp[— i(fj,s+eV)t/h], whereas the hole wave 
function oscillates as iph ~ exp[— i(fig — eV)t/%]. The 
noise combines these dependences in the product V'eV'h 
which now oscillates at the Josephson frequency 2eV/% 
corresponding to the singularity. In a junction containing 
a single superconductor, the singularity therefore appears 
in the noise rather than the current, and the detection 
of this frequency can be considered as an analog to the 
Josephson effect. 



B. Small biases: Andreev regime 

When the applied bias is much smaller than the gap, 
the S matrix elements can be taken to be constant. Us- 
ing the unitarity of the S matrix, both expressions (|l~6| ) 
and arc unified in the same formula over the whole 
energy interval < hu> < 2eV: 

4e 2 

S(uj) = —(2eV-hoj)R A (l-R A ) if htu < 2eV , 



. 5(w) = 



if Tiuj > 2eV , 
(18) 



where Ra = \sNNhe(Q)\ 2 is the Andreev reflection proba- 
bility. The noise spectral density decreases linearly with 
frequency, and vanishes beyond the Josephson frequency 
2eV/h (figure [l]), thus displaying a singularity at this fre- 
quency. 

This result has to be compared with both the Joseph- 
son effect £3 and with the analog result for a normal metal 
junctionEI (figure [l]). In the former DC bias ap- 

plied to a junction between two superconductors gen- 
erates an oscillatory current. The order parameter on 
each side oscillates as -01,2 ~ exp[— i2fj,g t 2 t/h] with fig 1 
and /xg, the chemical potentials of each superconductor 
(figure ^|a) . The resulting current involves the overlap 
of these two states V'iV'Ij an d therefore oscillates at the 
frequency 2\fig 2 ~ MsJ/ft- The noise characteristic ex- 
hibits a peak at 2&V/h which radiation line- width was 
computed in RefEa (inset figure |l]). In the former case 
(figure |b), the wave functions have a time dependence 
as -01,2 ~ ex p[ — iftiflt/h], so that although the resulting 
current is constant, finite frequency noise involves the 
overlap ipiipZ leading to a singularity at the frequency 
|A*2 - Mi I = eV/h. 

In the N-S case (figure §c), only Andreev reflection 
contributes to the current, involving the emission or the 
absorption of Cooper pairs (charge 2e) on the supercon- 
ducting side. An incoming electron in the normal side at 
energy [ig + eV drags another electron at energy ^g — eV, 
forming a reflected hole of energy /j,g — eV. The two 
combined electrons have a total energy 2fj,g which cor- 
responds to a Cooper pair, and are thus allowed to be 
transfcred to the superconducting side. Now the above 
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2eWh 




hco 

FIG. 1. Noise as a function of frequency. Full line: 
N-S junction, with a singularity at Tiui — 2eV (in units of 
(4e 2 / h)eV Ra(1 — Ra))', dashed line: junction between two 
normal metals with a singularity at Tiui = eV (in units of 
(4e 2 /h)eVT{l - T)). Inset: noise in the Josephson effect, 
with a peak at the frequency Tiuo = 2eV. 
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FIG. 2. Three different kinds of junctions: a) two su- 
perconductors in contact (Josephson effect), with chemical 
potentials /j,s 1 and fj,s 2 ; b) two normal metals with chemical 
potentials fii and fi2\ c) a normal metal ((j,n) connected to a 
superconductor (fis)- 



C. Large biases 

The applied bias is now larger than (or comparable to) 
the gap, so that the scattering matrix elements depends 
on the energy. A specific description is needed to charac- 
terize this dependence: the BTK modeled is particularly 



5 



suited for this purpose as it allows a description with a 
minimal number of parameters. 

A local tunnel barrier Vb(x) — VbS(x) is introduced 
at the boundary. Since the complete knowledge of the S 
matrix of the junction is necessary to compute the noise, 
the quasi-particle states in the superconductor are speci- 
fied. This only makes sense if these states are not evanes- 
cent (E > A). The Bogolubov-de Gennes equations for 
these states are solved on the superconducting side. Us- 
ing the continuity of the wave functions at the interface, 
and specifying the discontinuity of their derivative, one 
obtains the S matrix elements (Appendix ^). 

The energy integrals in Eqs. ( |l~6| ) and (17]) are per- 
formed numerically. Plotting the noise as a function of 
frequency, additional cusps or singularities are found at 
lu = (eV - A)/h, uj = (2A)/h, to = {eV + A)/h, on 
top of the Josephson singularity at u = 2eV/h (figure^). 
All these frequencies can be illustrated on an energy di- 
agram (figure ||) . This numerical calculation can also be 
performed for small biases, yielding full agreement with 
the previous calculation (Eq). Another interesting limit 
arises when eV 3> A. In this case, transport is dominated 
by single quasi-particle transfer with a charge e, whereas 
the contribution of Andreev reflection is small. Thus 
similar results to those of normal-normal metals junc- 
tion are expected. This is obviously the case (figure ||), 
even though the above mentioned singularities can still be 
identified. One may object that if the applied bias is too 
large, the non-equilibrium processes dominate, and the 
previous assumptions are not correct anymore because 
of heating effects. This limit is then valid for supercon- 
ductors with a small gap (A/ks ~ O.IK) because the 
condition eV A may be satisfied in a near-equilibrium 
situation. 

In order to visualize the additional singularities, the ar- 
gument invoking the oscillatory time dependence can be 
used once again. This time, because of the large value of 
the bias, several charge transfer processes occur: a) An- 
dreev reflection is still there, and it implies the same sin- 
gularity at the Josephson frequency 2eV/h. b) Electrons 
in the normal side are transmitted as electron-like quasi- 
particles in the superconductor. Wave functions oscillate 
as ipN,e ~ exp[— i(ns + eV)t/h] and ips,e ~ exp[— «(/ig + 
A)t/h], and the overlap gives a singularity at (eV — A) /Ti. 
Note that the same transfer process occurs with holes and 
hole-like quasi-particles, giving the same singularity, c) 
Electrons in the normal side are transmitted as hole-like 
quasi-particles in the superconductor (Andreev transmis- 
sion). Here, the time dependence of the wave functions 
is ip N , e ~ exp[-z(^ s + eV)t/h] and ips,h ~ exp[-i(^ s - 
A)t/h], implying a singularity at (eV + A)/Ti. The same 
transfer process exists with holes and electron-like quasi- 
particles. d) Andreev reflection also occurs on the super- 
conducting side, when electron-like quasi-particles are re- 
flected as hole-like quasi-particles, and vice versa. Wave 
functions oscillates as ijjs, e ~ exp[— i(ns + A)t/h] and 
ips,h ~ exp[— z(/i5 — A)t/h], giving a singularity at the 
frequency 2A/H. These four singularities are summarized 



in Fig. ft 




FIG. 3. Noise in a N-S junction as a function of fre- 
quency, with intermediate barrier transparency (Z = 1), for 
several values of the applied bias : eV = 0.8A, eV = 1.3A, 
eV = 1.9A. 
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FIG. 4. Energy diagram when the bias is larger 
than the gap. Relevant intervals associated with the 
cusps/singularities are shown. 
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where s ee , s e h. She and Shh are the standard matrix el- 
ements describing the N-S boundary only (without the 
external time perturbation). 
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FIG. 5. Noise in a N-S junction as a function of frequency, 
with intermediate barrier transparency (Z = 1), for a large 
value of the bias (eV = 20 A). 



D. Non-stationary Aharonov-Bohm effect 

Finite frequency measurements can represent a real 
challenge for noise, so it is interesting to imagine a sce- 
nario where an alternating field superposed to the DC 
bias allows to probe the finite frequency effect. The non- 
stationary Aharonov-Bohm effect has been introduced 
several years ago in a normal conductor connected to 
reservoirsB. In this proposal, a time dependent vector 
potential is applied in a confined region [xi,X2] °f the 
conductor, which adds a phase to the electrons and holes. 
The phase is chosen to be a periodic function of time 
<f>(t) = $ a sin(Ot) with $ a = 2ir dxA x /4> Ql and where 
4>o = he/ e is the normal flux quantum. The most striking 
consequence is the presence of steps in the derivative of 
the shot noise with respect to the voltage dSJ deV when 
the applied bias eV is a multiple of the frequency of the 
perturbation Hfl. Moreover, the gaps between the steps 
are non-monotonic with the amplitude of the harmonic 
vector potential. This effect has been experimentally ob- 
served in normal diffusive samplesEl 

Here, this result is extended to an N-S junction using 
the same framework (figure ||). The perturbation remains 
confined in the interval [xi, X2] near the boundary and is 
assumed to contain an adiabatic time modulation. The 
main difference with the previous case is due to the An- 
drecv reflection. The wave function of an incoming elec- 
tron accumulates a phase $(i) in the region where the 
potential is confined, but after a normal reflection the 
same phase is subtracted. On the contrary, if the elec- 
tron is Andreev reflected, the outgoing hole accumulates 
a phase totalizing a phase 2$ for the complete reflec- 
tion process. So in the Andreev regime, the S matrix of 
the junction can be written as follow: 



S = 



She e 



2i*(t) 



Seh e 

Shh 



2i$(t) 



+3> 



e <r 



+<f> 



h <^ 



(19) 



FIG. 6. The NS boundary. The light shaded region deter- 
mines where the electron/hole wave function may accumulate 
phase: schematic description of the two scattering processes, 
Andreev and normal reflection. 



In contrast to the usual Aharonov-Bohm effect, no 
closed topology is imposed: the flux does not have to 
be enclosed in a loop, and the current is not periodic in 
the accumulated phase 2$(i). Moreover, the effect of the 
perturbation on the average current is straightforward in 
the limit where the probability of Andreev reflection Ra 
depends weakly on the energy: it brings a periodic mod- 
ulation of the current AI = (4e 2 /h)RA[hH,/e]^ a cos(fit). 
The most striking consequence of the flux occurs when 
the current-current correlations are considered: the mod- 
ulation leads to a non-monotonic effect as a function of 
phase, in contrast with the electromotive force action on 
the current. 

However, because of this periodic perturbation, trans- 
lational invariance in time is broken, and the process be- 
comes non-stationary. So the noise depends in general 
on two frequencies and thus can be written as: 

S(Qi,Q, 2 )= [ I ' dt 1( it 2 e l{nitl+n2t2) 

{I(h)I(t 2 )) </) 2 ) . (20) 

This double Fourier transform can be rewritten as: 

+00 

S(fii,fi 2 )= 27rS ( n i + n 2 ~mn)s {m} (n 2 ) . (21) 

m— — 00 

The zero harmonic, which is proportional to S(fli + Q 2 ) 
is the most standard quantity to study. From Eq. Jl^), 
its zero-frequency expression is given by: 
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S<°>(0) 



2*2 



e 2 h 



1 



2m 2 v 2 F (2tt^) 2 „ ' —OO J — OO 



dt' 



+ OO /"+CO 

< J dE J dE lSr - 



f m {E){l-UE'))e^ E '- E ^'^' h 

(E,E',f) 

+E* NmNn {E, E ,t)B Nm N n (E, E', t) 

+ANmNn{E, E , t)B-NmNn{E, E , t) 
+ B NmNn( E 1 E ' ') 0^iVmJVn( £ i *') 

+/ m (S)/ n ( J B , )e- l(B+£ ' )(t '- t)/?l 

X C* NmNn (E, E') (CNnNm(E' , £7, t') 

(E,E' ,t') ) 
+ (1 - / m (E))(l - /„(£'))e l(£+£ ' )(t '^ )/?l 

X \ CN n N m (E' , E, t) + CNmNn{E, E' , t) 



XC 



n(E,E',t') 



(22) 



Performing this integral at finite temperature from the 
explicit time dependence of the current matrix elements 
and using the generating function of the Bessel functions 
J„, one obtains: 



+00 

£ 

rn— — 00 



4-p 2 

SW(p) = -j r R A (l-R A ) J] J^ a )F v (mhn) 



8r 2 

+ ^-R\k B T . 



(23) 



where the temperature dependence appears in the form: 

F v (mhtt) = (2eV - mhQ) cothpe^ - mhtt)/2k B T] . 

(24) 

Note the factor 2 reminiscent of the Cooper pair charge 
in the argument of the Bessel function, which originates 
from the accumulated phase 2<fi. Now the derivative of 
the noise with respect to the voltage is taken: 



+M 



" - ^R A (l-,l A ) ^ 4,2<T»„,. (25, 

m——M 



dV 



Fv(TnMl) specifies how the steps in the noise deriva- 
tive are smeared with temperature. In Eq. (|2^) the 
sum over hapnonics has a cutoff at M = [2eV/hil\ . In 
experiments)!!] it is more convenient to characterize the 
non-monotonic dependence on voltage by taking the sec- 
ond derivative of the Aharonov-Bohm contribution to the 
noise. This is illustrated for two distinct temperatures 
in figure ^[ For small temperatures k B T < hi}/ 2, the 



noise steps are individually resolved, and one observes 
oscillations as a function of 2eV/hCl, with a clustering 
of peaks with a large amplitude. For larger tempera- 
tures (k B T > Ml/2), one expects the signal to vanish 
but this is obviously not the case: although individual 
peaks separated by Kl can no longer be identified, clus- 
ters of peaks (or clusters of "large" steps in the noise 
derivative) continue to give an average contribution to 
the non-stationary AB effect. This robustness enhances 
the likelihood of experimental observations. 

The corresponding experiment waSp-successfully re- 
cently achieved in diffusive conductorsE-J For compar- 
ison with theory, the current spectral density is averaged 
over the transmission channels, and Eq. ( f23[ ) becomes: 

S(°)(0)=4fc B TG NS (l-v) 

+2 V G NS J2 J 2 m (^a)F v (mhn) , (26) 

m— — 00 

where Gns is the differential conductance of the device, 
and 77 = 1/3 the suppression factor for a normal diffusive 
conductor. The second derivative with respect to the 
voltage of the measured noise clearly shows peaks at the 
Josephson frequency. This constitutes a rather robust 
experimental check of the presence of an effective charge 
2e in the fluctuation spectrum of single NS junctions. 




0.0 2.0 4.0 



FIG. 7. Non-stationary Aharonov-Bohm effect: plot 
of d 2 S/dV 2 , expressed in units of (8e 4 /nh 2 n)R A (l - Ra), 
as a function of 2eV/hQ, with the choice $ a = 3. For 
2k B T = 0.27ifi (full line) and for 2k B T = HQ (dashed line). 



IV. HANBURY-BROWN AND TWISS 
GEDANKEN EXPERIMENT WITH A 
SUPERCONDUCTOR 

A. Introduction 

So far, attention on effective charges have been the pri- 
mary focus. Effects associated with the statistics of the 
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charge carriers are now considered. In the mid-1950s, 
Hanbury-Brown and TwissB described a new type of in- 
terferometer in order to find the size of a radio star by 
measuring the correlations between the signals of two 
aerials. This experiment was followed by another one 
using a coherent light source, where a mercury arc lamp 
beam was partitioned by a splitter into a reflected and 
a transmitted part. Intensity correlations between re- 
flected and transmitted beams were found to be positive. 
This result can be explained by the quantum statisti- 
cal properties of photons, which are bosons. Particles in 
a beam of bosons tend to cluster together (bunching), 
so the probability to detect simultaneously two photons 
(one in each beam) is non-zero, and therefore correlations 
are positive. On the other hand, two indistinguishable 
fermions exclude each other because of the Pauli principle 
(anti-bunching), and consequentki- reflected and trans- 
mitted beams are anti-correlatecffll. 

More recently, two analogs to the original Hanbury- 
Brown and Twiss experiment, using fermions propagat- 
ing in semiconductors (electrons) - instead of photons 
is vacuum - were achieved. NegatiV|S|-correlations were 
expected, and experimentally verified™! Here an hybrid 
system is envisioned: a superconductor is introduced just 
next to the beam splitter. This implies that transport in- 
volves Cooper pairs on the superconducting side, which 
have a finite penetration length, an the normal side, the 
essence of the proximity effectE3~E3. While Cooper pairs 
are not strictly bosons, an arbitrary number of these can 
exist in the same momentum state, so bosonic statistics 
could be detected in such a system. Thus the possibility 
for positive noise correlations cannot be ruled out. Be- 
low, it is shown that the sign of the noise correlations 
depends on the transparency of the beam splitter. 



5i 2 (0) 



2c 



2 r eV 



rev 

/ dE [ E w 



iee S ljeh s lihe S ljhh) 



s 2jeh s 2iee ~ S 2 jhh s ^ihe) 

+ (s liee Si4 ea — S li ^ e Si4^ Q ) 

i=l,2 ; a—e,h 

X ( s 24ea s 2iee ~ S^ha^ihe) ■ (28) 

The sign of Eq. (^8|) cannot be determined at this stage: 
it depends on the specific form of the S matrix. Thus, for 
analytical purposes, a simple model is proposed, where 
the beam splitter is dissociated from the N-S junction 
(see figure |g). To establish the S matrix of the junction, 
one needs to combine the S matrices of the beam splitter 
and of the N-S junction (between 3 and 4). 




FIG. 8. Two normal leads (1 and 2) are linked by a 
semi-transparent beam splitter (BS) and connected to a su- 
perconductor (4) via a normal region (3). 



B. Model 



C. S matrix of the splitter 



The device is composed of two normal leads (1 and 2, 
see figure ^]) linked by a semi-transparent beam splitter 
(BS) and connected to a superconductor. The state cor- 
responding to an electron incoming (outgoing) in (from) 
the lead i is labelled cf e (c£" e ). The hole incoming (out- 



going) in (from) the lead i is described by c i 



(see 



figure |3|) . With this convention, the S matrix of the whole 
system is defined as: 



t<he\ 

"2/i 



-4e 



s 



S4A 

c ih 

< 2 c. 

C 2h 
r .+ 
L 4e 

V C 4h ) 



(27) 



The goal is now to calculate the expression of the noise 
correlations between leads 1 and 2. From Eq. (15) and 
in the limit of zero temperature, it is possible to show 
that the expression for the noise correlations reduces to: 



The beam splitter is described by a S matrix Sm which 
gives the outgoing states as a function of incoming ones. 
For electrons one obtains: 



=2e 



"3e 



s 



M. 



■ -le 

n+ 

■ -2e 
,+ 
"3e 



(29) 



Expression of Sm c is identical to that of Ref.t 

a b y/e 
Sm. = | b a s/e 

y/e y/e -(a + b) 



(30) 



where a = (yl - 2e - l) /2, b = (y/l 
e can vary from to 1/2. Sm, 



2e + 1) /2, and 
depends only on a sin- 
gle parameter e which monitors the transparency of the 
splitter. For example if e = no transmission occurs 
from region (1) or (2) to region (3). A similar relation 
holds for holes: 



9 









( ^ } 






C 2h 


- Sm h 










\4J 




\ C 3h J 



(31) 



Note that the splitter does not couple electrons and holes. 
Expression of the matrix for holes is given by the relation 
SM h {E) = Sm (-E), as no magnetic field is assumed to 
be present here, but since Sm c is real and does not de- 
pend on the energy Sm h = Sm c ■ 



D. Small biases: Andreev regime 

When the applied bias is much smaller than the gap 
(eV <C A), Andreev reflection between 3 and_.4 is the 
only transmission process. One can then writeEj: 



<-3e 
c 3h 



7 
7 



>-3e 
-3ft 



(32) 



with 7 = e 4arccos ( B / A ). i n such a case, all matrix ele- 
ments like supq or S24 P q (with p,q — e, h) are zero. Set- 
ting x = VI — 2e it is possible to show that: 



Sllee — Sllhh — S22ee — S22hh 1 



S21ee — S21hh — S\2ee — S\2hh — " 



2(1 - 7 2 a; 2 ) 
(.T + l)(l-7 2 J r) 

2^2\ 



2(1 - 7 2 a; 2 ) 



(33) 
(34) 



Slleh — S21eh : 



Si2eh — S22eh — Sllhe 1 

7(1 -x)(l +x) 

2(1 -J 2 X 2 ) ' 



: S21he 



-Si2he — S22he 

(35) 



Because E <C A one can make the assumption that 
7 ~ —i. Since the S matrix elements do not depend on 
the energy anymore, the integral ( |2§| ) can be performed. 
One finally obtains: 



2e 2 e 1 



2e + 1 



(36) 



The noise correlations vanish at e = 0, when conductors 1 
and 2 are equivalent to a two-terminal device decoupled 
from the superconductor, and in addition, S12 vanishes 
when £ = \/2 — 1. A plot of S12 (normalized to the noise 
in 1 (or 2) at e = 1/2) as a function of the beam splitter 
transmission (figure indicates that indeed, the corre- 
lations are positive (bosonic) for < e < \[2 — 1 and 
negative (fermionic) for \[2 — 1 < e < 1/2. At maximal 
transmission into the normal leads (e = 1/2), the cor- 
relations give the negative minimal value: electrons and 
holes do not interfere and propagate independently into 
the normal terminals. This is the signature of a purely 
fermionic system. When the transmission e is decreased, 
Cooper pairs may leak in region 3 because of multiple 
Andreev processes at the boundary. Further reducing 
the beam splitter transmission allows to balance the con- 
tribution of Cooper pairs with that of normal particles. 



Expression (|36j) predicts maximal (positive) correlations 
at e = 1/3: a compromise between a high density of 
Cooper pairs and weak transmission. 
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FIG. 9. Noise correlations between the two normal reser- 
voirs of the device, as a function of transmission of the beam 
splitter, showing both positive and negative correlations. 



E. Larger biases 



If the applied bias is greater than the gap, transmission 
of quasi-particles between 3 and 4 is now allowed, and one 
has to take into account the en ergy d ependence of the S 
matrix elements. As in section [II C, the BTK model is 



chosen, and thus quasi-particles at arbitrary energy can 
be handled. For simplicity, the same beam splitter is 
still used, assuming that its S matrix is independent of 
the energy. Even though this may appear as a crude ap- 
proximation, this remains correct for example when the 
superconductor has a small gap. The S matrix of the 
whole system is computed in appendix |n|. 

A numerical calculation of noise correlations is per- 
formed with the help of Eq. (28). If one considers a high 
transparency barrier (Z = 0.1) and a small bias, one finds 
agood agreement with previous analytical results (figure 
|l0| ), except that the noise correlations (still normalized 
to the noise in a normal lead when e = 1/2) do not quite 
reach the minimal value at e = 1/2. This is an early sig- 
nature of the potential barrier at the NS interface. The 
more the bias is increased, the weaker are the positive 
correlations. If the voltage is large enough (beyond the 
gap) positive correlations are completely destroyed. As 
in the analogy with the normal-normal junction encoun- 
tered previously, Cooper pairs contribute only for a small 
part to the transport, and the system loses its bosonic 
features. 
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FIG. 10. Noise correlations using an N-S boundary mod- 
eled by BTK for a high transparency barrier (Z = 0.1). From 
top to bottom, eV/A = 0.5, 0.95, 1.2, 1.8. 

Intermediate transparencies are now considered {Z = 
1, figure [ll]), and a strikingly different behavior is ob- 
tained. For weak biases, noise correlations remain pos- 
itive over the whole range of e. It is possible to find 
an appropriate value (for example eV = 0.95A) in or- 
der to observe oscillations between positive and negative 
values of the correlations. Further increasing bias, corre- 
lations again become negative over the whole range of s. 
Calculations for larger values of Z confirm the tendency 
of the system towards dominant positive correlations at 
low biases with Si 2 (e)/5i(l/2) > lover a wide range ofe 
(not shown) . The phenomenon of positive correlations in 
fcrmionic systems with a superconducting injector is thus 
enhanced by the barrier opacity at the N-S boundary. 
Nevertheless, at the same time, for opaque barriers, the 
absolute magnitude of Si and S12 becomes rather small, 
which limits the possibility of an experimental check in 
this regime. 
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FIG. 11. Noise correlations using an N-S boundary mod- 
eled by BTK for intermediate transparency: Z = 1 (same 
biases as in Fig. hoj) . 



A suggestion for the experimental device is depicted 
figure |l2|. Assume that a high mobility two dimen- 
sional electron pgas has a rather clean interface with a 
superconductor^. A first point contact (Pi) close to the 
interface selects a maximally occupied electron channel. 
The beam of electrons is incident on a semi-transparent 
mirror similar to the oae,used in the Hanbury— Brown and 
Twiss fermion analoga§tl A second point contact located 
in front of the mirror (P2) allows to modulate the reflec- 
tion of the splitter in order to monitor both bosonic and 
fcrmionic noise correlations. In addition, by choosing a 
superconductor with a relatively small gap, one could ob- 
serve the dependence of the correlations on the voltage 
bias without encountering heating effects in the normal 
metal. 




^N 2 



FIG. 12. Proposed device for the observation of posi- 
tive/negative correlations; at the boundary of a superconduc- 
tor (5), two point contacts (Pi) and (P2) are connected to a 
semi-transparent mirror (M). 



An alternative interpretation of these positive correla- 
tion has been proposedcSE^ in terms of the sign of the 
effective charges obtained from the spectral weight as- 
sociated to electrons and holes. Nevertheless, statistical 
analogies are often useful in condensed matter physics as 
these allow to isolate the dominant behavior in the trans- 
port characteristic of a given system. For instance, in 
the fractional quantum Hall effect, dissipationless trans- 
port occurs because electrons with aa-pdd number of flux 
quanta represent a composite bosonEZI. 



V. CONCLUSION 

Both dynamical and statistical aspects of noise have 
been presented in a unified formalism. Expressions for 
the finite-frequency noise and the noise correlations for a 
conductor containing an arbitrary number of terminals, 
connected to a superconductor (Eq. ( fl4| ) and fll5|)) have 
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been derived. In a first step, a single N-S junction was 
considered. In the Andreev regime, the noise spectral 
density presents a singularity at the Josephson frequency 
and vanishes beyond. This can be interpreted in terms of 
an effective charge 2e transfered at the boundary. Note 
that this argument fails if the bias voltage is increased 
above the gap. In this case, both a single and a double 
charge transfer are allowed. Thus the effective charges 
which are identified in the spectral noise density plots 
are not well defined (e < e* < 2e). 

The non-stationary Aharonov-Bohm effect has been 
proposed as a tool to analyze these features in a zero fre- 
quency measurement. It allows the observation of peaks 
in the second derivative of the noise with respect to the 
voltage when the frequency of the applied perturbation 
is commensurate with the Josephson frequency. Because 
these kinds of measurements are in principle easier to 
achieve, this non-stationary Aharonov-Bjohm effect has 
caught the interest of experimentalistsEj'tj, who have 
provided a confirmation of the theoretical predictions of 
this paper. Thus, a single N-S junction is an adequate 
system to observe the Josephson frequency with only one 
superconductor instead of two, as in the usual Josephson 
effect. 

In a second step, the feasibility of an analog to the 
fermionic Hanbury-Brown and Twiss experiment with 
a superconductor has been addressed. Correlations are 
shown to be either positive or negative, depending on 
the reflection coefficient of the beam splitter. Therefore, 
such a fermionic system can exhibit a bosonic behav- 
ior. A qualitative interpretation of this puzzling result is 
reached as follows: when the transmission at the inter- 
face decreases, Cooper pairs leak on the normal side, and 
these can be considered as "composite bosons" , hence the 
positive statistical signature. Recent experiments in the 
"normal" fermionic Hanbury-Brown and Twiss analogpu 
were performed successfully, and these experiments could 
possibly be extended to the case presented here, giving 
the opportunity to observe for the first time positive cor- 
relations in a fermionic system. 

Nevertheless, the issues of this paper present some lim- 
itations. All the calculations have been performed in 
the single channel case, even if a generalization to sev- 
eral channels is possible. However, this extension should 
not bring major changes. The assumption that the S 
matrix elements are independent on the energy is cor- 
rect as long as the applied bias remains small enough. 
For larger biases, the BTK model takes this energy de- 
pendence into account, but only allows numerical calcu- 
lations. Moreover, further increasing the bias leads to 
a non-equilibrium situation which cannot be described 
with this formalism. An appropriate approach would be 
to employ the Keldysh Green's functions method. How- 
ever, in order to go further than the present calcula- 
tions, heating effects should be taken into account self- 
consistcntly. Nevertheless, the physics presented here is 
rather robust, and different results, especially concern- 
ing the effective charge 2e in the Andreev regime are not 



expected. The calculations have been performed for a 
arbitrary S'-matrix .describing a specific sample, but as 
pointed out in Ref.tLZl, all these results can be extended 
to the diffusive case without difficulty by averaging over 
the ttansmission channels with the standard methods of 
RefB. 

Future considerations may include the transport char- 
acteristics of other types of superconductors (high T c 
superconductors), where the gap varies in momentum 
space. Recently, the noise in a junction between a 
normal metal and a d-wave superconductor has been 
calculatecO. Moreover, as it was emphasized above, 
the study of the statistics and of the effective charges 
of quasi-particles is a relevant issue in condensed matter 
systems. For instance, in the fractional quantum Hall 
effect, Laughlin's ewasi-particles are supposed to obey 
fractional statistical. Results concerning particles which 
obey exclusion statisticscil showed that the shotpthermal 
noise crossover deviates from the fermion casec3. An- 
other example is the study of the transport properties 
of "atomic" composite fermions/bosons, such as alkali 
atoms or 3 He and 4 He which off-equilibrium properties 
are not known at this time. Finally, noise correlations 
have been computed here only at zero frequency. Is there 
anything to gain from the finite frequency spectrum of 
noise correlations ? This quantity combines the dynam- 
ical aspects of current fluctuations with the statistical 
nature of multi-terminal geometries. 
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APPENDIX A: S MATRIX ELEMENTS OF N-S 
JUNCTION WITH A BARRIER 

Following the BTK modelH, with the barrier potential 
Vb(x) = VbS(x) and the pair potential A(x) = A6(x), 
the Bogolubov-de Gennes equations are solved on each 
side of the junction, giving the states in the normal lead 
and in the superconductor. The wave functions are con- 
tinuous at the N-S boundary. Integrating the Bogolubov- 
de Gennes equations on each side of the junction gives 
another condition over the derivatives of the wave func- 
tions. Writing these conditions for both particle types 
and both sides of the junction, one obtains four linear 
systems of equations with four unknowns, which are the 
S matrix elements: 

{ul-vt){Z* +l Z) 

SNNee = SsSee = (Alj 

7 



12 



SNNhh = SsShh 



{ul-vl){Z*-iZ) 



SNNhe = SNNeh = —SSShe = — SSSeh 



UqVq 

7 



SSNee = SNSee 



SSNhh — SNShh 



u oV u o- v o( 1 + iZ ) 



(A2) 
(A3) 

(A4) 

(A5) 



SSNhe = —SSNeh = —SNShe = SNSeh 



where: 



1/ (E 2 -A 2 )V 2 



= 2 (1 



= 1 



(A6) 



(A7) 



Z = tjiVb I (h 2 kp) is the relative height of the barrier and 
7 = Ug + (mq — Vq)Z 2 . Note that the unitarity of the S 
matrix is satisfied. 



APPENDIX II: S MATRIX ELEMENTS OF A THREE TERMINAL N-S JUNCTION WITH A BARRIER 



To compute the complete S matrix, one searches which outgoing states are obtained when a single particle is 
injected in a given terminal. For example, injecting an electron in 4 (ie a state c^ e ), one obtains reflected waves (c^ e 
and c£ h ) and transmitted waves (cj~ e , c 



matrix elements: s 4 4 ee , suhe s 14ee , s 14he 
types of particle. One obtains: 



c 2e and c+ 



S24e 



(see Fig. |8|). One can therefore deduct the corresponding S 
and S2ihe- This operation is made for all the terminals and for both 



S44ee S44 e ^ 
S44he Stehh 



~>SS 



(a + b)S S N[l+ (a + b)S NN ] 1 S NS , 



(1) 



s 14ee Si4 e /j 



S24ee S24eh 
S24he S24hh 



~(S N s - (a + b)S NN [1 + (a + b)S NN ] 1 S NS ^j 



(2) 



S41ee S4i e h 
S41he S4ihh 



S42ee S42e/i 
S42he S42hh 



= VeSsN [1 + (a + b)S NN ] 



(3) 



Sllee 
Sllhe 



S22ee 
S22he 



a 




e + a(a + b) 
e 



S 



NN 



[1 + (a + b)S NN ] 



(4) 



Sllee 
Sllhe 



S22ee 
S22he 




a 



£ 

e + a(a + b) 



S. 



NN 



[l + ia + ^SNN]- 1 f 



(5) 



Sl2ee 8i2eh 
Sl2he Si2hh 



S21ee S21eh 
S21he S21hh 



Sllee + 1 Slle/t 
Sll/te Sllh/l + 1 



(6) 



Here Sns, Ssn, Snn and Sss are 2x2 matrices of which elements have been computed in appendix [X|. All the 
calculations have been made above the gap. But the obtained results are valid even if the energy is smaller than the 
gap, using the correct values for the S matrix of the N-S junction: Sns and Ssn are zero, and Snn and Sss are 
given in appendix |A|. In this case, expressions of Sn, S22, S21 et S12 remain the same, but 541, S42, S14 and S24 are 
zero, and S42 = Sss- 
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